Reorientation of the large-scale circulation in turbulent Rayleigh-Benard convection 
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We present measurements of the orientation 9o(t) of the large-scale circulation (LSC) of turbulent 
Rayleigh-Benard convection in cylindrical cells of aspect ratio 1. 9o(t) undergoes irregular reori- 
entations. It contains reorientation events by rotation through angles A8 with a monotonically 
decreasing probability distribution p(A8), and by cessations (where the LSC stops temporarily) 
with a uniform p(A0). Reorientations have Poissonian statistics in time. The amplitude of the LSC 
and the magnitude of the azimuthal rotation rate have a negative correlation. 

PACS numbers: 47.27.-i, 05.65.+b, 47.27.Te 



A major component of turbulent Rayleigh-Benard con- 
vection 0, 0, E| (RBC) in a fluid heated from below 
is the emission of hot (cold) volumes of fluid known 
as "plumes" from bottom (top) thermal boundary lay- 
ers. These plumes are carried by, and by virtue of their 
buoyancy in turn drive, a large-scale circulation (LSC) 
[3- 0- IB- HI known as the "mean wind" . For cylindri- 
cal samples of aspect ratio T = D/L ~ 1 (D is the sample 
diameter and L the height) the LSC is a circulation in a 
near-vertical plane with up-flow and down-flow on oppo- 
site sides close to the side wall. An interesting aspect of 
the LSC observed by others is that it can spontaneously 
and erratically change directions 0, 0, 0, 0] , either 
by an azimuthal rotation of the entire structure with- 
out much change of the flow speed, or by reversing its 
flow direction without much rotation of the circulation 
plane. We shall refer to the latter process as cessation 
since it involves a momentary vanishing of the flow speed, 
and to both processes collectively as reorientations. Ro- 
tation through an angle AO = 1/2 (we measure AO in 
units of 2ir, i.e. of revolutions) was observed in exper- 
iments but to our knowledge cessation so far had 
been demonstrated only in numerical simulations of a 
two-dimensional system llOi and in a dynamical-systems 
model of the phenomenon [14j. Thus, the re-orientation 
of the LSC dynamics is still only partly understood. 
Aside from its fundamental interest, the phenomenon is 
important for instance because it occurs in natural con- 
vection of the Earth's atmosphere 0] , and because it is 
responsible for changes in the orientation of the Earth's 
magnetic field that is a result of convection reversal in 
the outer core [lfi |. 

In this letter we describe experiments that shed new 
light on the nature of the reorientation of the LSC. We 
found that there is a broad spectrum of events leading 
to reorientation. At one extreme there is relatively slow 
rotation, over a time period of order ten turnover times 
T of the LSC, without a significant decrease of the circu- 
lation speed. At the other, there are near-instantaneous 
cessations that involve a vanishing of the flow speed fol- 
lowed by a re-initiation of the circulation with different 
angular orientations 0q. The re-orientations by either 



mechanism correspond to a continuous range of AO, in- 
cluding the rotation by AO ~ 1/2 observed before [ll] ] 
and the cessation involving AO ~ 1/2 that was seen in 
two-dimensional simulations lQj (where only this is pos- 
sible), favored in the interpretation of some experiments 
[l^j l , and contained in a recently developed model 0] . As 
seen by Sreenivasan et al. [hlj l , the temporal succession of 
reorientations had a Poissonian distribution, suggesting 
that successive events are independent of each other. 

The experiments cover the range 3xl0 8 ^i?5;10 11 of 
the Rayleigh number R = agATL 3 /nv (a is the isobaric 
thermal expansion coefficient, g the acceleration of grav- 
ity, AT the applied temperature difference, k the thermal 
diffusivity, and v the kinematic viscosity). The Prandtl 
number was a = v / 'k = 4.4. Two cylindrical samples 
with aspect ratio T ~ 1 were used and are the medium 
and large sample described in detail elsewhere ^3] • Both 
had circular copper top and bottom plates and a plexiglas 
side wall. The medium (large) sample had D = 24.81 
(49.67) cm and L = 24.76 (50.61) cm. Each was filled 
with water at a mean temperature T m = 40.0°C. Eight 
blind holes, equally spaced azimuthally in the horizontal 
mid-plane of the samples, were drilled from the outside 
into the side walls. Thermistors were placed into them 
so as to be within 0.07 cm of the fluid surface. They de- 
tected the location of the upflow (downflow) of the LSC 
by indicating a relatively high (low) temperature, but did 
not perturb the flow structure in the fluid. We made mea- 
surements with a sampling period St ~ 6 seconds, and 
fit the function = T + S cos(i7r/4 - ), i = 0, . . . , 7, 
separately at each time step, to the eight temperature 
readings. Here S is a measure of the amplitude of the 
LSC and 0q is the angular orientation of the plane of cir- 
culation. Typically the uncertainties were about 13% for 
5 and 0.02 for . 

In Fig. we show a long time series of 0q from the 
large cell and for R ~ 10 11 . It covers a period of nine 
days. In this example, on average one sees net rotation, 
at a rate of about one revolution per day in the clock- 
wise direction when viewed from above. Although 6o(t) 
samples all values of 0, there is an orientation m where 
the probability distribution p(0q) has a maximum. Net 
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FIG. 1: Time series of the orientation 8o (in revolutions) and 
of the amplitude S (in °C). (a): long time series of Go from the 
large sample for R — 9.4 x 10 10 . (b) to (d): time series during 
reorientations from the medium sample for R = 1.1 ■ 10 10 for 
(b): a rotation; (c): a reversal; and (d): a half-reversal. 



rotation was observed also by Sun et al. |18| . but in the 
counter-clockwise direction. 

Hidden away in time series like that in Fig. ^ (a) are 
numerous reorientations. A few examples for the medium 
cell and R = 1.1 ■ 10 10 are shown in Figs. ^ to d (for 
this case T = 49s) . Figure Q (b) shows a net angular 
change AO ~ 1/2 over a period of about 10T, similar 
to an event reported in Ref. [ll|. The reversal time 
much larger than T indicates a relatively slow rotation 
of the LSC orientation. Throughout this event 6 is seen 
to remain non-zero. 

Figure ^ (c) also shows an angular change of approx- 
imately half a revolution, but the change occurred in a 
fraction of T. Here S decreased for about one turnover 
time before the reversal, reached approximately zero dur- 
ing the reversal, and then increased for about another T. 
One sees that the mean wind gradually slowed to a stop, 
reversed direction, and gradually sped up again without 
significant rotation of the plane of circulation. This ces- 
sation is the mechanism favored by Sreenivasan et. al. 
|l3| in the interpretation of their data. Presumably it 
would prevail in systems where the rotational invariance 
is broken, either by imperfections of a nominally cylindri- 
cal sample or by a deliberate other choice of the geometry 
such as a square horizontal cross section. As we shall see 



below, in our system such events are relatively rare. 

In Fig. (d) the plots are qualitatively similar to 
(c). As in (c), the vanishing of 5 indicates that the LSC 
stopped (also a cessation), but when it re-started, it as- 
sumed an orientation that differed by A9 ~ 1/4 from the 
original one. Many other values of A9 were observed, 
indicating that there is a larger class of phenomena with 
a continuous distribution of AO. 

In order to automatically identify members of the ex- 
pected large spectrum of reorientation events, we re- 
quired reorientations to satisfy two criteria. First, the 
magnitude of the net angular change in orientation \A9\ 
over a set of successive data points for 8 had to be 
greater than A9 m i n . Second, the magnitude of the net 
average azimuthal rotation rate |(0o)| = \A9/At\ over 
that set had to be greater than 9 m m- Here At is the 
duration of the reorientation. Often there were mul- 
tiple overlapping sets that satisfied these requirements. 
Then the set with the maximum local quality factor 
Q n = |At9|/(At) ra was chosen as the reorientation. For 
< n < 1, Q n represents a compromise between choos- 
ing the maximum angular change (Qo) or the maximum 
rotation rate (Qi). Any adjacent points to the chosen 
set were also included if the rotation rate <9 = 59o/St 
[S9 = 9o(t + St) — 9o(t)] for the additional range was 
greater than 9 m i n and of the same sign as the reorien- 
tation. Mostly we used the parameters A9 min — 0.125, 
dmin = 0.1 /T, and n = 0.25; but since all three are 
arbitrary, we did the analysis over the ranges 1/64 < 
A9 min < 0.25, 0.0125 < 6 min T < 0.2, and < n < 1 to 
confirm that the results are not sensitive to the choice. 
For A9 m in and 9 m %n , the smallest values used were close 
to the noise level of the data, and above the largest values 
too few reorientations were counted to yield useful statis- 
tics. The qualitative conclusions of the analysis did not 
change with different parameter values. However, some 
of the quantitiave values changed. This will be mentioned 
as it is appropriate. 

An obvious feature of reorientations is the angular 
change A(9. All of the reorientations found for the 
medium sample were sorted into bins according to At9|, 
regardless of the Rayleigh number R. The probability 
distribution p( | A9 1) is shown in Fig. |2|with statistical er- 
ror bars, in which the relative error is equal to the inverse 
square root of the number of reorientations in each bin. 
Fitting all data for p(\A9\) by the maximum-likelihood 
mehod to a power law in the range |At9| > 0.125 
yielded p(\A9\) ex (|A6»|) 7 , with 7 = -3.74 ± 0.04. The 
same analysis was done also for reorientations at various 
fixed R. These probability distributions also could be 
fitted by power laws, with the resulting 7 values shown 
for each R in the inset of Fig. |3 One sees that the 
distribution was, within our resolution, independent of 
R. The power-law exponents were approximately con- 
stant over large ranges of the reorientation definition pa- 
rameters; however they could vary by a factor of 2 for 
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FIG. 2: Probability distribution of the angular change p(A8) 
for reorientations. Solid circles: All data. Solid line: power- 
law fit to the data. Inset: power-law exponent 7 versus R 
for the medium cell (solid diamonds) and large cell (open 
diamonds). Dashed line: the value of 7 from the fit to all 
data. Note that reversals with A9 ~ 0.5 are relatively rare 
with p(A8 = 1/2) ~ 0.1, implying that only about 1 % of the 
events correspond e.g. to A9 = 1/2 ± 0.05. 
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FIG. 4: Solid circles: probability distribution p{t\/(t\)) of 
the time intervals n between succesive reorientations, (a) over 
the entire range of Ti measured, and (b) for ri < (ti). Open 
circles: the probability distribution p(Ari/(ri)) of the dif- 
ference between successive intervals. Solid lines: the func- 
tion p(ti/(ti)) = exp(— ri/(n)) representing the Poisso- 
nian distribution for p(n/(ri)). Dashed line: the function 
p(Ati/(ti)) = 0.5exp(— |Ati/(ti)|) representing the Poisso- 
nian distribution for p(Ari/(ri)). Dotted line: fit of a power 
law with exponent —1 to the data. 



FIG. 3: Probability distribution of the angular change p(\ A8\) 
for cessations. Here |A#| is reduced to the range < |A#| < 
0.5. Solid circles: All reorientations with (5) < 5(i?)/4. Solid 
line: the uniform distribution. 



extreme values of the parameters. Qualitatively in all 
cases, the distribution was consistent with a power law 
for \A9\ > 0.125 with a large negative i?-independent ex- 
ponent. The major qualitative conclusions are that there 
is a continuous distribution of AO, that smaller reorienta- 
tions are much more common than larger ones, and that 
there is no characteristic reorientation size. 

We would like to distinguish between cessations and 
rotations, so the probability distribution of the angular 
change AO for cessations only is plotted in Fig. [31 There 
is not always a clear distinction between cessations and 
rotations, but here we used the criterion that a reori- 
entation is a cessation when the average value (S) dur- 
ing the reorientation satisfied (S) < 8{R)/A, where S(R) 
is the average amplitude of an entire run at a given R. 
This probability distribution for cessations is consistent 
with a uniform distribution, which is very different from 
the distribution for all reorientations shown in Fig. [5] 
Since cessations make up only a small percentage of all 
reorientations, the probability distribution for all events 
as shown in Fig. |3 still accurately describes rotations. 
The uniform distribution of angular changes for cessa- 



tions means that at the end of a cessation the orientation 
of the LSC is equally likely to be at any angle. This sug- 
gests that the mean wind loses its memory of its initial 
orientation during the cessation. 

Another important aspect of reorientations is their dis- 
tribution in time. All of the time intervals t± between the 
ith and (i + l)th reorientations from the medium sample 
were sorted into bins according to the value of T\. The 
probability distribution p(ri/ (ti)) is shown in Fig. El(a), 
with statistical error bars as before. The average time in- 
terval between reorientations (ti), computed separatly at 
each R, was approximately 307", but quantitatively de- 
pended on the definition of reorientations and on R. The 
data are in good agreement with the exponential func- 
tion p(ti/{ti)) = exp(— ti/(ti)), which represents the 
Poissonian distribution. 

The Poissonian nature of reorientations can be seen 
further by considering the difference between successive 
time intervals. Thus, Ari = Ti i+i — was sorted into 
bins according to At\/ (t±) and the probability distribu- 
tion p(A"n/ (ti)) is also plotted in Fig. 0] (a). The data 
are in good agreement with the function p(Ati/(ti)) = 
0.5exp(— |Ati/(ti)|), which again represents the Poisso- 
nian distribution. 

The analysis in terms of t\ and Ati implies that succes- 
sive reorientation events are independent of each other. 
For large T\ / (n ) this is consistent with the conclusion 
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FIG. 5: (a): The negative cross-correlation g^ s (r) between the 
magnitude of the rotation rate \0\ and the amplitude 5 at R — 
1.1 • 10 10 (solid circles) and a fit by an exponential function 
(solid line).(b): Reynolds numbers for the average r 9i o [solid 
(open) circles] and width r Sl i [solid (open) diamonds] of gg S (r) 
for the medium(large) sample. 



of Sreenivasan et. al. for reversals in their experi- 
ment that were defined very differently. However, these 
authors also found that for small time intervals, with 
Ti < 30T in their experiment, that p(r{) followed a power 
law distribution with an exponent of — 1 . To better com- 
pare the present work with this result, p(t\/{ti)) is plot- 
ted in the range t\ < (n) ~ 30T in Fig. |U(b). Our data 
in this range are also consistent with the Poissonian dis- 
tribution as found for the larger range, and inconsistent 
with the —1 power law found by Sreenivasan et. al. This 
difference between the two experiments is likely due to 
the different methods of defining and analyzing events. 
We tested this by alternatively defining a reorientation 
as any crossing by 9o(t) of an angle orthogonal to the pre- 
ferred angle 9 m of the LSC orientation. This definition is 
more similar to that of Ref. 01 • I* 1 that case we obtain a 
p(ri/ (ti)) that rises well above the Poisson distribution 
as n/ (ti) decreases below one and that can be described 
quite well by a power law with an exponent close to one 
for ri/(n) < 1. 

Finally we consider the relationship between the mag- 
nitude of the azimuthal rotation rate \6q\ and amplitude 
6, since the examples in Fig. 2] show \9q\ is large when 
S is small. A negative cross-correlation function between 
the two quantities is given by 

9et(r) = -((\e (t)\-(\d \))(S(t-T)-S(R))}/(a § a s ); 

#o(i)|-(N)) 2 >; 
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In Fig. (a) we show gg S (r) for R 
The function gg S (r) = a exp(— |(t 



1.1 • 10 
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? ,o)/r g ,i|) was fit 
to the data separately at each R to yield r gj0 (-R) and 
T 9t i(R). These time scales were made dimensionless 
by dividing the viscous diffusion time L 2 /v by each of 
them, yielding Reynolds numbers i?f'° = L 2 /(r g fiv) and 
R9^ = L 2 /(T g sv). These are plotted versus R in Fig. 



□ (b). A power-law fit gives Rf° = 0.41 • R 0A1 and 
Rf 1 = 0.051 • i? 46 . One sees that these two power 
laws have approximately the same exponent. This expo- 
nent value is found also for the Reynolds number for T 
|2fj| , suggesting a relationship between the different time 
scales of the problem. The time scale t 9i \ indicates how 
long the rotation rate and amplitude remain correlated, 
and it is about 80% of T. The time scale r g fl indicates 
by how much the amplitude leads the rotation rate. This 
is only about 6% of T, but the fact that the amplitude 
leads the rotation at all is remarkable, and suggests that 
an amplitude drop is a precursor of reorientations. 
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